Topological jamming of extended structures and the glass transition 
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We propose a new scenario for glassy dynamics in frustrated systems with no quenched-in ran- 
domness, based on jamming of extended dynamical structures near a critical point. This route to 
a glassy state is demonstrated in a lattice model of fluctuating lines. Numerical simulations of the 
model show non-exponential relaxations and diverging energy barriers in the vicinity of a thermo- 
dynamic phase transition. A master equation for the coarse grained dynamics is constructed. It 
shows how topological jamming leads to the observed glassy dynamics. 



Foams, granular materials and supercooled liquids, ex- 
hibit a transition from a flowing liquid phase to a jammed 
phase 0. The dynamics of the flowing state near the 
transition are characterized by slow relaxations mediated 
by large-scale spatial rearrangements. These dynamical 
heterogeneities were dramatically demonstrated in recent 
experiments on colloids near the glass transition An 
important question posed by these observations is: What 
is the relation between length and time scales in a glass 
forming material? 

The mode-coupling theory predicts a divergence of re- 
laxation time scales with no dramatic growth of length 
scales. This is born out by experiments. A scaling theory 
of the glass transition, based on the existence of a zero- 
temperature critical point leads to a similar conclu- 
sion about the relation between length and time scales. 
The frustration-limited domain theory relates the ori- 
gin of the dynamical heterogeneities and the activated 
nature of the dynamics to a critical point "avoided" due 
to frustrations imposed by competing interactions. Sim- 
ilarly, the random first-order transition theory ||^ relies 
on the presence of frustration in glass formers and sug- 
gest that an unusual critical point is responsible for the 
dynamical behavior observed near the glass transition. 

A phenomenological jamming phase diagram was re- 
cently proposed as a means of putting under one umbrella 
the whole spectrum of systems exhibiting dynamical ar- 
rest 1^. How, and if the above mentioned theories relate 
to jamming ideas remains to be explored. Here we pro- 
pose a scenario for glassiness based on the interplay of 
jamming and an underlying critical point. We demon- 
strate it in detail in a lattice model of fluctuating lines. 

a. Loop model: In constructing a lattice model aimed 
at describing supercooled liquids, we concentrated on two 
key features; frustration and the presence of extended 
structures. A model which has these attributes is the 
three coloring representation of the loop model |^]. The 
model is defined on a 2-dimensional hexagonal lattice, 
in which all bonds are colored A, B or C with the con- 
straint that no two bonds which share a vertex are of the 
same color. The number of colorings is exponential in 
the number of vertices and the system has an extensive 
entropy. This entropy is maximized when the number of 
bonds of any color, in each of the three directions is the 
same. This entropy was calculated exactly by Baxter 
A coloring with all the bonds in one direction having one 



and the same color has zero entropy density. By intro- 
ducing an energy parameter that favors colorings in the 
zero entropy sector a phase transition can be induced in 
the model. 

The frustration of the coloring model gives rise to nat- 
urally occurring extended structures in the system. Con- 
sider all the A and B colored bonds: there is one of 
each at every vertex, and therefore they must form loops. 
All the properties of the model, including the correla- 
tions of the microscopic variables, the colors, can there- 
fore be traced back to the behavior of the ensemble of 
loops. These loops do not intersect, they have orienta- 
tion {ABAB... versus BAB A...), and therefore they 
can be thought of as contour lines of a height field, that 
lives on the dual lattice and provides an appropriate 
order parameter for the loop model. 

The lattice with periodic boundary conditions, has a 
linear dimension L with 2L^ sites and 3L^ bonds and the 
three lattice directions are indexed by a = 1,2,3. We 
introduce a long-range interaction energy between C col- 
ored bonds, with a parameter measuring the strength 
of this interaction: 
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Here is the total number of C bonds in the a 



di- 



rection. In the ground state all of the C colored bonds 
are aligned in one of the three directions. This ordering 
effect acts against the entropic influence which tries to 
keep the number of colors homogeneous in all directions, 
and hence disordered. As the value of the coupling con- 
stant /i is increased from 0, the system undergoes a phase 
transition from a disordered to an ordered phase. 

The phase transition is best described in terms of the 
height representation. In the disordered state, where the 
colors are randomly arranged, the corresponding surface 
is rough and untilted; Fig. |l|(a) shows the height of a typ- 
ical coloring state. In the ordered state, with no exten- 
sive entropy, the surface is smooth and maximally tilted 
as shown in Fig. 0(6). 
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FIG. 1. The height profiles in the (a) disordered, and (b) 
ordered states are shown, for L = 24; the x— and y— axes 
denote the lattice coordinates, and 2— axis gives the height. 



example, a ...ACAC... loop becomes a ...CACA... loop, 
with probability 1 if the energy change is non-positive 
and with rate exp(— Ai?) when Ai? is positive. Since the 
energy function defined in Eq. ^ depends on global num- 
bers and not on local ones, it is easy to see that all local 
loop updates keep AE = 0. Only when system span- 
ning loops with nonzero winding numbers are chosen for 
update, does the possibility of A£' ^ arise. There- 
fore, the tilting of the surface and hence the lowering of 
energy can only be brought about by updates of system- 
spanning loops. This path downhill in energy is blocked 
by barriers which lead to spectacularly slow dynamics 
near /x*, as will be shown below. 



c. Simulation results: We measure the equilib- 
rium, macroscopic tilt-tilt autocorrelation function, cor- 
responding to the slowest mode in the system: 
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A tilting transition of the surface occurs at /i = 
The existence of a transition at a finite value, /z* can be 
argued as follows. The entropy is maximum for a rough, 
untilted surface |^ . Similarly from Eq. |l| it follows that 
the energy is also maximized at zero tilt. Therefore the 
free energy, F — E — S , which is concave at fi — 0, must 
change its curvature at a finite /i» . We numerically mea- 
sured the equilibrium probability distribution of tilts for 
different values of ^, and observed that the distribution 
becomes broader with increasing /i. Quadratic fits indi- 
cate a /i* ~ 0.35 ± 0.05, at which the second derivative 
of the free energy function at zero-tilt vanishes. A tilting 
transition at a finite fi^ has been explicitly demonstrated 
in the two-color version of the current model 1^ lending 
further support to our argument. 

b. Loop dynamics: The tilting transition, whose dy- 
namics we are interested in studying, involves a change 
in the density of non-zero winding number loops: struc- 
tures which span the system. The existence of these 
mesoscopic structures and the role they play in the tilt- 
ing transition is reminiscent of the extended dynamical 
structures observed in atomistic simulations |10| and ex- 
periments near the glass transition. This motivates 
us to consider loop dynamics. 

In a unit microstep in our simulation, we randomly 
choose a site and choose either an A — B , B — C 01 a, C — A 
loop passing through it, and we switch the colors on ev- 
ery alternate bond, using ordinary metropolis rules. For 



Here p = {px,Py) is the tilt vector, and (• • •) denotes 
an average over various initial times to in a history. We 
obtain an average C(t) by averaging several ('^ 10) histo- 
ries. The time is measured in units of Monte-Carlo step 
per site (MCS). In Fig. ^ we show C{t) at different p for 
L — 24. Similar curves have been obtained for two other 
system sizes, L ~ 36 and L = 48. The autocorrelations 
functions exhibit a two step-decay, which may be termed 
the (3— and a— relaxations following the standard termi- 
nology of mode-coupling theory | ll| . Although we do not 
attempt a detailed analysis of the shape of these relax- 
ation curves in this letter, several features are apparent 
from Fig. ^ with increasing /i the time scale of relaxation 
to the plateau, the plateau height, and the time scale of 
the decay of the plateau increase monotonically. In fact 
the decay time scale of the plateau r (see inset of Fig. 
H), grows exponentially and diverges at a finite /z with 
the divergence being well described by the Vogel-Fulcher 
form; exp(A(L)/(^*(L) — fi)) for all three system sizes 
studied. The extrapolation of fJ,*{L) to the infinite sys- 
tem indicates a divergence close to the transition point 
/i, obtained from the equilibrium distribution. 
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FIG. 2. The log-log plot of C{t) versus t, At fi = 0.14, 0.15, 
0.175, 0.2, and 0.22 for L = 24. The curves exhibit a 2-step 
decay, and the fits (shown as solid lines) provide us with an 
estimate of the decay time scale of the plateaus, r. The inset 
shows a plot of r against ^ and a Vogel-Fulcher fit to it. 

As discussed in detail in the following paragraphs, the 
exponential divergence of t and the shape of C{t) can 
be associated with jamming barriers and enhanced static 
fluctuations near the critical point. 

d. Residence time and transition rates of different 
tilted states: The Vogel-Fulcher divergence of the relax- 
ation time is a generic feature of real glass formers. The 
origin of this anomalous slowing down is still a matter of 
some debate. Within our simple model we investigate the 
crucial features underlying the rapid rise in time scales, 
by measuring the transition probabilities among different 
tilted states and constructing a coarse-grained dynamical 
model for the tilt relaxation. We find that the relaxation 
process is very different from simple Langevin dynamics 
within a concave free energy well and is, instead, domi- 
nated by traps and barriers. 

Although the autocorrelation functions C{t) were mea- 
sured with boundary conditions which allow the system 
to have any tilt p, to simplify our discussion, we concen- 
trate on simple tilted states with the tilt pointing either 
in the x or y direction. The measurements primarily in- 
volve studying several histories (of the order of 10"* — 10^) 
in which the system starting from an initial state p, stays 
in that state for time t and makes a transition to a new 
state p in the next time interval St. From these time 
histories, one can easily calculate Pp^p, the probability 
that the system stayed in the state p for time t + St, and 



Pp~tp' J the probability that it stayed in the state p for 
time t and made a transition to p between t and t + St. 

Assuming that the tilt configurations {p} evolve via 
a Markov process, the measured probabilities should 
have an exponential form: Pp^p = /iexp(— t/rp) and 
Pp^p' = /2 exp(— t/rp), where /i and /2 are numbers 
less than unity. From our extensive simulations we find 
that at each p and for any p (with p ^ 0) the probabil- 
ities Pp^p and Pp^p' (for say p = {p + 1) and (p — 1)) 
have a common exponential factor with identical decay 
constants Tp. This evidence validates the above descrip- 
tion in terms of the coarse-grained states {p}. The p — 
states are an exception since we found that Po^o is a 
stretched exponential in t. This observation has little 
bearing on our current discussion, and will be taken up 
in the future. 

Based on the above arguments, the measurements of 
Pp^p and Pp^p' , can be used directly to calculate the 
relaxation times Tp and the transition matrix elements 

Wp'p from p to p : 

p p p 

pVp 



The variation of the relaxation time scales with the tilt 
p are shown in Fig. ||(a) for different values of p. The 
most striking feature of these time scales is that Tp in- 
creases exponentially with p: Tp oc exp(Q;(p^''')) with 7 
between 1/2 and 1. We will show below that the expo- 
nential increase of Tp arises from the jamming of non-zero 
winding number loops and is directly responsible for the 
Vogel-Fulcher divergence of the average time scale. 

To explore the origin of the exponential rise of Tp with 
p, we investigated the elements of the transition matrix 
W obtained from our measurements. The rates for the 
transition Wp_i,p, and Wp+i^p are plotted in Fig. ^(6) 
and Fig. ^(c), respectively for different p's as a function 
of p. 




(b) (0) 



3 



FIG. 3. We show (a) the exponential dependence of Tp on 
p for fj. — 0.15, 0.175, 0.2 and 0.22. The major contribution 
to the plateau timescale r in C{t) {cf Fig. ^) comes from 
non-simple tilts. For example, Tp for p = (4,2) is ~ 211 and 
660 MCS, for p = 0.15 and 0.175, respectively. (b) IVp-i.p 
versus /x at p = 1, 2, 3 and 4 (from top to bottom), and (c) 
Wp+i,p versus /j. sd, p — 1, 2 and 3 (from top to bottom). 

Transition rates from higher to lower tih are well fitted 
by exp (— 8/z(2p — 1)). This suggests these rates to be of 
the form: 



W, 
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(4) 



The fj, independent factor To varies rather weakly (about 
a factor of 2 for the data in Fig. ^(6)) with p, and thus 
the relaxation of p to {p — 1) is dominated by the cost of 
going uphill in energy as expected. 

The unusual behavior of Tp can be traced to the tran- 
sition rates Wp+i^p (Fig. ||(c)) which can be written as: 

W^p+i,p = F<,F(p). (5) 

The "blocking factor", T{p), shows a negligible depen- 
dence on p but a strong, exponential dependence on the 
tilt p. The origin of T(p) can be traced back to jamming 
caused by the system spanning loops in a tilted state. 
These loops are meandering objects, and their lengths 
scale as L^/^ [ p"2|Jl3| thus squeezing in an extra loop in 
order to increase the tilt is associated with straightening 
out the existing ones. Thus the transition to higher tilt 
states are dramatically suppressed with increasing p and 
lead to an exponential increase of Tp. Using Eqs. ^, ^ 
and 1^, the time scales can be expressed as: 
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The p-dcpendent blocking factor in the denominator of 
Eq. cannot be scaled away and the dramatic suppres- 
sion of the the downhill transitions, reflected in T{p) <C 
exp(— (i?p_i — Ep)) < 1 leads to a Tp being dominated by 
the exponential energy cost of decreasing the tilt 

e. Theory of Vogel-Fulcher divergence: We have es- 
tablished that the relaxation times Tp grow exponentially 
with p because of barriers arising from the spatial con- 
straints imposed by the system-spanning loops. In this 
paragraph we discuss how the shape of C(t) and the 
Vogel-Fulcher type divergence of r can be related to the 
Tp. We have good numerical evidence that the time de- 
pendence of C{t) is well approximated by the form: 
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This form of C{t) is based on the approximation that 
the system is perfectly correlated as long as it stays in a 
state with a given tilt and decorrelates completely as it 
leaves this state. The probability of staying in a tilt p for 
a time t is measured by e"*/'^" and Peq{p) measures the 



equilibrium probability of the tilt p. The above model of 
C{t) is identical to that of the trap model of glasses [p^ . 
Since Tp grows exponentially with the tilt ~ exp(ap^'''), 
the plateau timescale in C{t) is dominated by the Tp cor- 
responding to the maximum accessible tilt, Pmax, as de- 
termined by the width of the distribution Peq{p). As the 
transition to the tilted state is approached, the fluctu- 
ations of p increase leading to a power-law divergence 
of (yO^) and Pmax and an exponential divergence (Vogel- 
Fulcher) of the plateau timescale, t. 

In summary, we flnd that a lattice model with frustra- 
tion and associated extended structures exhibits a trap 
dominated dynamics arising due to jamming. A phase 
transition involving these structures then naturally leads 
to an exponential divergence of timescales. 
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